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ABSTRACT 



Context. The determination of lithium abundances in stars of young clusters have shown that they deplete Li by different degrees 
during their pre-main sequence phase. Beryllium abundances are complementary to the lithium ones, and can help tracing the mixing 
processes in the stellar interiors. 

Aims. Our aim is to derive beryllium abundances in a sample of G- and K-type stars of two young pre-main sequence open clusters, 
IC 2391 and IC 2602. The Be abundances are used to investigate the mixing of internal material in these stars. The reliability of the 
Be lines as abundance indicators in low-temperatures is also investigated in detail. 

Methods. We derived Be abundances from high-resolution, high signal-to-noise UVES/VLT spectra using spectrum synthesis and 
model atmospheres. Atmospheric parameters and other elemental abundances are adopted from a previous work. 
Results. The sample stars have masses in the range between 0.80 < M/M Q < 1.20. They have been shown to differ in lithium 
abundance by about 0.60 dex, with lower A(Li) in cooler and lower mass stars. Here, we find that all the stars have the same Be 
abundance within the uncertainties. These observations show that the Be abundance is not affected by the mixing events in the pre- 
main sequence, in this mass range, in agreement with the expectation of evolutionary models. A comparison with Be abundances in 
older clusters shows that, contrary to the models, cool stars deplete Be during their main-sequence lifetime, confirming what has been 
previously suggested in the literature. 
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1. Introduction 

Open clusters are ideal targets to study stellar evolution and stel- 
lar structure. They are composed of stars with the same age, 
chemical composition, and distance; moreover, their masses can 
be relatively well estimated. In young clusters, these advantages 
can be used to investigate the pre-main sequence (PMS) evolu- 
tion of low-mass stars. 

In this context, the evolution of the light elements lithium 
and beryllium offer valuable information. These elements burn 
in (p,a) reactions at relatively low but different temperatures in 
the stellar interiors (~ 2.5 x 10 6 K for Li and ~ 3.5 x 10 6 K for 
Be). Variations in their surface abundances, with respect to their 
initial values, help constraining the physical pro cesses mixing 
the surface material to that in the stellar interior |pTnsonneaultl 
ll997tlDelivannis et al]l2000tlCharbonnel et alJl200fj ~ 

Li abundances have, consequently, been extensively studied 
in open clusters. According to standard evolutionary modelfl Li 
depletion should be a unique function of age, mass, and chemical 
composition, in other words, within a given cluster stars with 
similar effective temperature should have undergone the same 
amount of Li depletion. Observations, however, have revealed 
a complex behavior with important discrepancies in almost all 



evolutionary stages of low- and intermediate-mass stars (see e.g. 
iPinsonneaul j| 1 997l and references therein). 

Focusing on PMS Li depletion, several surveys of young 
clusters have evidenced two major discrepancies between theo- 
retical predictions and empirical patterns. First, standard models 
seem to predict more Li depletion than observed in PMS clus- 
ters (defined as clusters with age between ~ 10 and 100 Myr). 
Second, there is a significant spread in Li abundances among 
late-G and K dwarfs, which has first been detected in t he 120 
Myr Pleiades dDuncan & Joneslll983trButler et al.lll987h . along 
with a connection between Li abundances and rotation; namely, 
fast-rotating stars appear to have higher abundances than slow- 
rotating ones. This connectio n was confirmed in a la rge sam- 
ple of Pleiades members by ISoderblom et alJ d 19931) and has 
been found to hold for stars betwee n 0.7 and 0.9 M Q , but to 
break down for less massive stars dGarcia Lopez et all 1 1994b 
I Jones et~al l ll996h . The scatter in Li abundances in otherwise 
similar stars has since then been confirmed by other authors 
(I.Teffrieslll999HFord et al.ll200llKing et al"1l2010h . 

Similar results of a star-to-star scatter in Li and the 
Li-rotation connection have been found in other clusters 
and associations both younger and older than the Pleiades, 
e.g. the nine young associations with ages between 6 and 



* Based on observations made with the ESO/VLT, at Paranal 
Observatory, under program 66.D-0284. 

1 By standard evolutionary models we mean those that only allow for 
mixing in convective layers and do not include effects due to rotation 
and/or magnetic fields. 



70 Myr analyzed in Ida Silva et al.l (12009b. IC 4 665 (~ 28 
Myr, Martin & Montesl 119971: I Jeffries et al.l l2009h . IC 2602 
(~ 45 Mvr. iRandich et al.l 1 19971 1200 ll) . NGC 2451 A and 



B (~ 50 Myr. iHunsch et al l 12004 a Per (~ 90 Myr, 



iBalachandran et alJll988l I20LU IRandich et alJll998l) . M 35 ( 
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175 Mvr.lBarrado v Navascues e"tafll2001ah. NG C 6475 (~ 220 
Myr, [Fames & Jeffries! 1 1997b ISestito et al.ll2003l) . and M 34 (~ 
250 Mvr. lJones et al.lll997l) . The general picture seems to indi- 
cate that the scatter decreas es with age and disappears by the age 
of the Hyades (~ 600 Mvr. lSoderblom et alj[l99a [19951) . 

The interpretation of the Li scatter as a real abundance 
scatter and its connect ion to rotation is sometimes disputed i n 
the literature (see e.g. iKing et alJ.l2000t IXiong & Dengj|2005t) . 
It has been argued that the Li scatter is only an apparent ef- 
fect caused by our failure of properly modeling the atmo- 
sphere of these stars, neglecting effects caused by stellar ac- 
tivity for example. Theoretical and observational investigations 
of these effects, however, seem to indicate that they are not 



TabFe 1. Physical data of the open clusters as adopted from the 
literature. 



able to fully explain the observed scatter ( Stuik et al. 1997t 



Barrado v Navascues et afl l2001b[ iRandichl 
20101) 



2001b IKing et al.l 



We finally mention that the observed trend of Li with ro- 
tation is the opposite of that predicted when effects induced by 
rotation, such as meridional circulation and shear tu rbulence, are 
included dChabover et al.ll995HMendes et al.ll999l) . In the mod- 
els, the fast-rotating stars are expected to deplete more Li then 
the slow-rotating ones. In the last 15 years several mechanisms 
related to rotation, magnetic fields, disk locking, and angular 
momentum loss have been proposed with the ai m of solving the 
discrepancies with observed PMS Li depletion ([Martin & Clare j 
1996b 1 Ventura et al. | ll998b IChabrier et all 120071: |Bouvierll2008b 
Baraffe & ChabrierlboiOt iMacDonald & Mullanll2010h~ WhiIst 



some of these models do predict a larger amount of Li depletion 
for slow rotators, they still need to be further investigated and 
constrained. 

Beryllium, on the other hand, has been investigated in far 
fewer stars than Li. This is because Be abundances can only be 
determined from the Be n 2 S- 2 Po resonance lines at 3131.065 A 
and 3130.420 A, a spectral region near the atmospheric cutoff 
and thus strongly affected by atmospheric extinction. Moreover, 
this near-UV region is extremely crowded with atomic and 
molecular lines, some of them still lacking proper identification. 

Most of the work o n Be in open cluster has concentrated in F 
and earlv-G stars (e.g. [Stephens et alj|1997b iBoesgaard & Kind 
2002b IBoesgaard et al.ll2003b[ 12004b iRandich et alJl2002[ 12007 



Smiljanic et al. 1 120101) . while very little has been done in cooler 
stars (r e ff < 5400 K). Stars with temperatures down to ~ 5000 
K have been investigated onl y in the Hyades (~ 6 00 Myr) by 
iGarcia Lopez et alj d 1995b andlRandich et al.l d2007l) and in NGC 
2516 (~ 150 Myr) bv lRandich et alJd2007l) . 

These works have shown the following behavior of the 
Be abundances. As in the case of the Hyades, although the 
Li abundance decreases by two orders of magnitude with de- 
creasing temperature bet ween 6000 to 5400 K, th e Be abun- 
dances remain unaltered dBoesgaard & Kindl2002l) For cooler 
stars in the Hyades, t he results of Garcia Lopez et all d 19951) and 
IRandich et all d2007l) suggest a small degree of Be depletion. For 
the four young stars in NGC 2516, no Be depletion was detected 
down to ~ 5100 K dRandich et al.ll2007l) . For older field stars, 
ISantos et alj d2004l) have found some degree of Be depletion in 
stars cooler than 5600 K. These studies therefore suggest a low 
level of main-sequence Be depletion in the cooler stars. 

Ex cept for one star of IC 2391 analyzed by IRandich et alJ 
d2002l) . Be has never been systematically investigated in cool 
late G- and K-type stars of young pre-main sequence clusters. 
Here, we make for the first time such an attempt, and extend 
the investigation of Be to stars with T e g as low as ~ 4700 K in 
the PMS clusters IC 2391 and IC 2602. As evolutionary mod- 
els fail to reproduce the Li depletion pattern in these clusters 



Cluster 


(m - M) 


E(B - V) 


Age 


Distance 


[Fe/H] 




(mag.) 


(mag.) 


(Myr) 


(pc) 




IC 2391 


5.80 


0.01 


50 


144.9 


-0.01 


IC 2602 


5.86 


0.04 


46 


148.6 


0.00 



Notes. The distance modulus and physical distanc es were determined 
using Hi pparcos parallaxes bvlvan Lee uwen (2009). The color excesses 
are from lPatten & SimorJ ill 9961) and lProsser et all d 19961) for IC 2391 
an d IC 2602, respectively. Th eir chemical composition was determined 
bv lD'Orazi & Randichl d2009h . The ages are discussed in the text. 



(IRandich et al.ll200"ll) . it is important to investigate Be in order 
to better constrain the properties of the mixing during the PMS 
phase. 

Our work is divided as follows. In Sec.|2]we describe the tar- 
gets and the observational data. In Sec. [3] we present the model 
atmosphere and abundance analysis, discussing in detail the dif- 
ficulties of modeling the spectra of rotating cool stars. Our find- 
ings are presented in comparison to evolutionary models and to 
literature data on older clusters in Sec. [4] Finally, in Sec. [5] we 
present our final remarks. 



2. Targets and observational data 

2.1. The young clusters IC 2391 and IC 2602 

IC 2391 and IC 2602 are among the closest c lusters to the 
Sun, with distances of the order of ~150 pc dvan Leeuwenl 
120091) . Although these two southern clusters are close to 
the Galactic plane, their proximity results in small red- 
denin g. The clusters also have very similar ages. For IC 
2391. iBarrado v Navascues et al.l d2004l) derived an age of 50 
Myr using the lith ium depletion boundary method (see e.g. 
iBildsten et aDll997l) . For IC 2602, the Li dep letion method has 
recently been applied by iDobbie et al.l d2010h and an age of 46 
Myr was found. The clusters are young enough that 1 M stars 
have just arrived or are arriving at the zero age main-sequence 
(ZAMSfl The data on the clusters is summarized in Tab. [TJ The 
young ages and proximity make the two clusters of considerable 
interest for the study of pre-main sequence evolutionary phases. 

In spite of that, only a few analyses of their chemical com- 
position h ave been publi s hed. T he metallicities were first deter- 
mined by IRandich et al.l d2001l) for a subsample of the early-G 
and late-K type stars that they analyzed. Values of [Fe/H]3 = 
-0.03 + 0.07 and [Fe/H] = -0.05 + 0.05 were f ound for IC 2391 
and IC 2602, respectively. IStiitz et al.l d2006l) analyzed 5 B-F 
stars in IC 2391 including star SHJM2, one of the stars included 
in our sample. A mean metallicity above solar of [Fe/H] = +0.06 
was determined, if we exclude HD 74044 which was found to be 
a slightly chemically peculi ar A-type sta r . For S HJM2 a value of 
[Fe/H] = +0.08 was found. iPlatais et al.l d2007l) derived an aver- 
age value of [Fe/H] = +0.06 ± 0.06 from high resolution spectra 
of four slowly-rotating G-dwarfs. More recently, detailed abun- 
dances for 15 G- and K-type dwarf members o f the two clus- 
ters were obtained bv lD'Orazi & Randichl d2009). Solar average 



2 According to the models presented in lPiau & Turck- Chiezd d2002h . 
for example, stars with 1.4 M G and 1.0 M Q , with Pleiades chemical 
composition, arrive at the ZAMS after pre-main sequence lifetimes of 
~30 and ~50 Myr, respectively. 

3 FA/B] = log [N(A)/N(B)]* - log [N(A)/N(B)] 
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metallicities of [Fe/H] = -0.01 ± 0.02 and [Fe/H] = 0.00 ± 0.01 
were determined for IC 2391 and IC 2602, respectively. Here we 
adopted these last values. 

2.2. Targets 

In this work we analyze six G- and K-type stars of each clus- 
ter (liste^Jn_Jabj2]^All_the stars were included in the work 
of iD'Orazi & Randichl (120091) . The color-magnitude diagrams 
of the clusters are shown in Fig. [1] where the observed stars 
are shown as open circles. The sample stars have been iden- 
tified as photometric counterparts of possible lo w-mass cluster 
members detected in ROSAT X-ray observations (iRandich et al.l 
119951; iPatten & Simon! 1 19961) . The membership of these stars 
is supported by spectroscopic obs ervations of radial veloci- 
ties, Ha activity, and Li abundances dRandich et al.Hl 99711200 it 



IStauffer et alJT997tlMermilliod et al.l2009l) For IC 239 1 , proper 
motions have been measured by Platais et afj d2007h . These mea- 
surements confirm the membership of all our sample stars, ex- 
cept for star VXR31, which is found to be a non-member. A 
search in the literature shows that this star has always been 
flagged as having a discrepant radial velocity, but was suspected 
to be a single-lined spectroscopic binary. We therefore assume 
here that it is not a member of the cluster, and we exclude it 
from further discussion. 

2.3. Observational data 

The observations were carried out in visitor mode using UVES, 
the U ltraviolet and Visual Echelle Spectrograph dDekker et alj 
|2000T) fed by UT2 of the VLT. UVES is a cross-dispersed echelle 
spectrograph able to obtain spectra from the atmospheric cut-off 
at 300 nm to ~ 1 100 nm. UVES was operated in dichroic mode 
with cross dispersers #1 and #4 and central wavelengths 346 nm 
and 750 nm, respectively. A slit of 0.8" was used, resulting in 
a resolution of R ~ 55 000 in the blue arm. The reduction was 
conducte d with the ESQ UV ES data reduction pipeline within 
MIDAS dBallester et al.ll2000h . The spectra have average signal- 
to-noise (S/N) between 55 and 115 in the region around the Be 
lines. The log book of the observations is given in Table|2] 



3. Analysis 

3.1. Atmospheric parameters 

The atmospheric paramet ers were adopted from the work of 
ID'Orazi & Randichl (I2009I) . Their analysis was done in LTE and 
differentially to the Sun. Effective temperatures were determined 
with the excitation equilibrium of Fe i lines. A value of log g = 
4.50 was in general adopted for the surface gravities, except for a 
few slow -rotating stars where the gravity could be checked using 
the ionization equilibrium between Fe i and Fe n lines. In these 
cases a final value of log g = 4.45 was found, in excellent agree- 
ment with the initial assumption. The microturbulence was con- 
strained by requiring the abundances calculated from Fe i lines to 
have a null correlation with the equivalent widths. The reader is 
referred to their work for more details on the adopted line list and 
atomic data. For the Sun, adopting the parameters, T e g = 5777 
K, log g = 4.44, £ = 1.10 km s" 1 , a value of A(Fe) = 7.52 was 
derived. As mentioned above, average metallicities of [Fe/H] = 
-0.01 ± 0.02 and [Fe/H] = 0.00 + 0.01 were determined forIC 
2391 and IC 26 02, respectively. The atmos pheric parameters are 
given in Tab. [3] ID'Orazi & Randichl (120091) estimate that the ran- 



dom errors affecting their parameters are of + 60 K for T e f[, + 
0.20 km s for and between + 0.15 and + 0.25 dex for log g. 

3.2. Be abundances 

The beryllium abundances were determined us ing synthetic 
spectr a calculated with the codes described in ICoelho et al.l 
(120051) . Model atmospheres were computed for the exact at- 
mospheric parameters of the stars us ing the Linux version 
(ISbordone et al.l 12004 ISbordond 120051) of the ATLA .S9 code 
originally developed by Kurucz (see e.g. lKuruc9l 1993b . For the 
ca lculations we adopted the new opacities distribution functions 
of ICastelli & Kuruczl (120031) without overshooting. These mod- 
els assume local thermodynamic equilibrium, plane-parallel ge- 
ometry, and hydrostatic equilibrium. 

The line list is the same as used in our previ ous works on Be 
abundances dSmiljanic et al. |[2008l[2009l l 2010h . The molecular 



line list is that described in ICoelho et al 



2005|^ and the atomic 



line list is the one compiled bv lPrimas et al.l (119971) . Log gf of 
-0.168 and -0.468 were adopted for the Be lines at 3131.066 A 
and 3130.421 A, respectively. The line list includes a Fe 1 line at 
3131.043 A, with log gf = -2.517 and X = 2.85 eV, with the ob- 
jective of simulating an unknown line that affects the blue wing 
of the Be 3131 line. The parameters of this line were constrained 
using several stars of different parameters and metallicities (see 
iPrimas et al.ll 19971 for details). The proper identification of this 
line is still controversial in the literature, and we remark that 
other choices hav e been made by differ ent authors. As discussed 
in more detail by IRandich et al] d20 07). this choice does not af- 
fect the calculation of Be abundances of stars with r e ff > 5400 
K. As a number of our stars are cooler than that we discuss this 
issue in detail below. 

With the line list described above, the solar UVES spec- 
trurrQ, and adopting the parameters: r e jf = 5777 K, log g = 4.44, 
and £ = 1 .00 km s _1 an abundance of A(Bejf| = 1 . 10 is obtained. 
This abundance is in exce llent agreement with the one found by 
IChmielewski et all d 19751) . A(Be) = 1.15, and usually adopted as 
the reference photospheric solar abundance. 

As for som e stars in the open cluster IC 465 1 analyzed in 
ISmiljanic et alj (1201 Oh . the spect ra of our sampl e stars are af- 
fected by rotational broadening. ISmilianic et al.l d2010l) found 
that from this kind of data reliable, albeit slightly overestimated, 
Be abundances can be determined. The reader is referred to the 
discussion presented in lSmilianic et al ] d2010h for more details. 
Examples of the fits for the stars VXR70 (from IC 2391) and 
R66 (from IC 2602) are shown in Figs. [2] and [3] respectively. 

We also derived Be abundan ces in two stars of the Hyades 
analyzed in IRandich et al.l (120071) . vbl09 and vbl82. The abun- 
dance we find for both of them is A(Be) = 0.95 + 0.23 (the 
two stars have temperatures of 5141 and 5079 K, respectively). 
That should b e comp ared with A(Be) = 0.80 + 0.21 derived in 
IRandich et all d2007D . 

3.3. The Be 31 31 line below 5400 K 

Of the two Be resonance lines, the one at 3130 A is stronger but 
also more affected by blends. Because of that, the line at 3131 
A is usually preferred for abundance determinations. The latter 
is, however, also affected by blends. The most important blends 



4 The spectrum is available for download at the ESO website: www. 
eso.org/observing/dfo/quality/UVES/pipeline/solar_spectrum.html 

5 A(Be) = log N(Be) = log(Be/H) + 12. 
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Table 2. Observational data. 



Star(l) 


Other 


V 


(B-V) 


(V-R) 


Date of 


Exp. Time 


S/N at 




Name 


(mag.) (2) 


(mag.) 


(mag.) 


observation 


(s) 


3130 A (3) 


VXR3a 


CD-52 2467 


10.95 




0.37 


16/02/2001 


2 x 2400 


100 


VXR31 


CI* IC 2391 PMM 6478 


11.22 




0.37 


17/02/2001 


3 x 2400 


20 


VXR67a 


CI* IC 2391 PMM 5859 


11.71 




0.52 


16/02/2001 


4 x 1800 


55 


VXR70 


CD-52 2524 


10.85 




0.38 


17/02/2001 


2 x 2400 


115 


V AK / L 


1 it ojoy-ZjU-1 


1 1 .40 






1 c/no nc\c\\ 
1 o/UZ/ZUUl 


4 X loUU 


/u 


SHJM2 


CPD-52 1614 


10.30 


0.57 


0.33 


16/02/2001 


2 x 1200 


90 


Rl 


TYC 8960-1957-1 


11.57 


0.91 




16/02/2001 


4 x 1800 


75 


R14 


TYC 8964-606-1 


11.57 


0.87 




17-18/02/2001 


4 x 1800 + 1 x 1650 


75 


R15 


TYC 8964-73-1 


11.75 


0.93 




18/02/2001 


5 x 1800 


60 


R66 


CPD-63 1624 


11.07 


0.68 




16/02/2001 


2 x 1800 


80 


R70 


HD 307926 


10.92 


0.69 




18/02/2001 


2 x 2400 


80 


R92 


HD 308013 


10.26 


0.67 




16/02/2001 


2 x 1200 


70 



Note s. The stars from IC 2391 are listed in the upper part of th e table and stars from IC 2602 in the lower part. (1) Identi fiers for IC 2391 are from 
iPatten & Simor] ill 9961 VXR#) and fromlStauffer et al.ldl989L SHJM#). For IC 2602 they are fr omlRandich et al.l dl995h . (2) The VR photometry 
of the IC 2391 members is from lPatten & Simorid 19961) and the (B-V) color of SHIM2 is from lStauffer et al.Hl989l) . The BV photometry of the 
IC 2602 members is from lProsser et al.l d 19961) and Randich et al. ( 1995). (3) Average S/N around 3130 A in the combined spectrum. 

1 1 1 1 I 1 1 1 1 I 1 1 
10 - IC 2391 

f ' 

>14 - ' "■.*■ • 



16 - 



18 - 

I I I I I I I I I I I L 

0.5 1 1.5 2 0.5 1 1.5 2 

(V-R) (B-V) 

Fig. 1. Color magnitude diagrams of IC 2391 (left panel) and I C2601 (right panel). The observed stars are shown as open circles. 
T he VR photometry of IC 2391 is from lPatten & Simoril dl996l) . Only membe rs and suspec t ed me mbe rs according to the an alysis 
of IPatten & Simonl d 19961) are shown. The BV photometry of IC 2602 is from lProsser et al.l d 19961) and lRandich et al.l (Il995l) . The 
photometry was obtained by means of the WEB DA database. In the CMD of IC 2391, star VXR31, found to be a non-member, is 
shown as a full triangle. 

are: (i) a Zn line at 3131.109 A at its red wing (included in 
our a nalysis with parameters determined by ICorliss & BozmarJ 
[19621) ; (ii) two OH (A 2 £-X 2 IT) lines and on e CH (C 2 S-X 2 IT) line 
(included with parameters determined by ICastilho et al.l 1 19991) 
at its blue wing; (Hi) an unk nown line at its blue wing that we 
model as a Fe i line following lPrimas et al.l dl997l) . 

The lower the effective temperature the weaker the Ben 
lines are expected to become, as the population of the ionized 
species decreases. The blending features from neutral atoms and 
molecules, however, increase in importance. One might then 
question down to what temperatures the Be lines are still reli- 
able abundance indicators. 

To address this issue, the behavior of the 3131 Be line in 
lat e-type stars cooler than 54 00 K was inv estigated in de t ail firs t 
by IGarcia Lopez et al.l d 19951) and then bv lRandich et al.l d2007l) . 
iGarcia Lopez et alj d 1995b adopted a different solution for the 
unknown blending line, i.e. increasing the log gf value of a Mn i 



_j_L 



10 



12 



14 



1 Q 




line at 3131.037 A by 1.50 dejfl This Mm line is also i ncluded 
in our line list, but with the orig inal log g f determined by Kuruczl 



d 19881) . i.e. log gf = -1.725. IGarcia Lopez et all d 19951) con- 
cluded that below ~ 5200 K, the Mm becomes the dominant 
component of the blend, therefore the utility of the 3131 Be line 
as a good abundance indicator decreases. The exact temperature 
where it might become useless, however, would depend on both 
the resolution and the S/N ratio of the spectrum. Because of this, 



6 We note that other tentative identificati ons for the blend h ave been 
advocated in the literature. For example, iKing et aTl d 19971) decided 
to enhance the log gf of the Mnn 3131.015 A line by +1.726 dex. 
However, as also pointed out in the paper, they failed to fit the spectrum 
of Procyon, since this choice predicts a feature that is too strong. In 
addition, from our exercise we conclude that a feature which increases 
in strength with decreasing temperature is more adequate. The ionized 
line of Mn would have the opposite behavior. 
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Table 3. Atmospheric parameters and abundances of Li and Be. 



Star 


7W (K) 


log g 




v sin/ 


A(Li) 


A(Be) 




(K) 




(km s" 1 ) 


(km s" 1 ) 






VXR3a 


5590 


4.45 


1.15 


10 


2.80 


1.30 


VXR67a 


4680 


4.50 


0.70 


8 


2.30 


1.30 


VXR70 


5557 


4.50 


1.15 


17 


2.80 


1.30 


VXR72 


5257 


4.50 


1.20 


15 


2.70 


1.10 


SHJM2 


5970 


4.45 


1.20 


< 15 


2.90 


1.20 


Rl 


5050 


4.50 


1.15 


< 10 


2.40 


1.25 


R14 


5100 


4.50 


1.10 


13 


2.65 


1.20 


R15 


4770 


4.50 


0.70 


7 


2.40 




R66 


5590 


4.45 


1.15 


12 


2.80 


1.30 


R70 


5760 


4.45 


1.10 


12 


2.90 


1.15 


R92 


5630 


4.45 


1.20 


14 


2.93 


1.20 


vB109 


5141 


4.50 


0.86 




0.49 


0.95 


vB182 


5079 


4.50 


0.84 




0.73 


0.95 



Notes. The atmospheric parameters for the stars of IC 2391 and IC 2602 we re derived bv lD'Orazi & R andich ( 20 0]|). In the bo t tom w e include two 
star from the Hyades with parameters determined bv lRandich et al. d2007t). The rotation al velocities are from lStauffer et all dl997l) . The lithium 
abundances for the stars of IC 2391 and IC 2 602 were calculated by Randich et al .1 J200 lh and include NLTE corrections. For the Hyades stars the 
Li abundances are from Randich et al. ( 2007). All the beryllium abundances were determined in this work. In addition, its spectrum has quite low 
S/N ratio and thus the value of A(Be) is quite uncertain. 





Fig. 2. Fit to the region of the Be lines in the star IC 2391 VXR 
70. The solid line represents the observed spectrum, the dashed 
line a synthetic one with A(Be) = 1.30, and the dotted lines spec- 
tra with ± 0.20 dex in the Be abundance. 



Fig. 3. Fit to the region of the Be lines in the star IC 2602 R 66. 
The solid line represents the observed spectrum, the dashed line 
a synthetic one with A(Be) = 1.30, and the dotted lines spectra 
with ± 0.20 dex in the Be abundance. 



iGarcia Lopez et ail (1 19951) were able to put only upper-limits to 
the Be abundance in their sample stars with T e ff < 5200 K. 

In lRandichetalJ d2007h . a similar exercise was conducted. 
These authors, however, adopted the same atomic line list that 
we use in this work and compared the di fferent effects caused 
by adopting either the Fe i or the Mn i line. lRandich et al.1(l2007l) 
concluded that for stars above 5400 K choosing either of the so- 
lutions does not affect in a significant way the Be abundances. 
For lower temperatures, they found that the Be abundances can 
differ by 0.15-0.10 dex between 5250-5400 K, and up to 0.20 
dex at ~ 5000 K. Smaller abundances are derived when the Fe i 
line is adopted (i.e., this line is foun d to be always s tronge r than 
the Mn one). Using the Fe blend, iRandich et al.l (120071) con- 
cluded that it was possible to measure Be down to the lowest 
r e ff in their sample, ~ 5000 K. 



Here, we repeat the exercise, comparing the Fe and Mn lines, 
but also investigating the influence of the Zr i and the molecular 
lines on the blend and extending the test down to 4800 K. A 
comparison of the relative strengths of these lines is shown for 
five different temperatures in Fig. [4] 

As expected, the Be n line becomes weaker with decreasing 
temperature while the lines of Fei, Mm, and Zn all become 
stronger. At ~ 5100 K, the added contribution of the blending 
lines are comparable to the strength of the Be li ne, irrespective 
of the choice between Fe or Mn. As was found bv lRandich et alJ 
(I2007h . the Fe line is always stronger than the Mn line. We 
also confirm the differences in the Be abundances caused by 
adopting different li nes in the temperatures range studied by 
IRandich et al.l (120071) . In addition, we find that for even lower 
temperatures, down to ~ 4800 K, the difference can amount to 
0.40-0.45 dex. 
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Fig. 4. Comparison of the strength of the lines affecting the Be line at 3131 A. The Be line is shown as a solid black line, the 
molecular features as a dotted black line, the Mn line (with increased log gf) as a green short-dashed line, the Fe line as a red 
long-dashed line, and the Zr line as a blue dot-dashed line. Five different temperatures, corresponding to the temperature range of 
our sample, are shown: (a) 5760 K, (b) 5590 K, (c) 5300 K, (d) 5100 K, and (e) 4800 K. 



While we have no definitive answer to whether the blend 
is caused by a Fe or a Mn line, it is clear that had we chosen 
to use the Mn line in our analysis, we would obtain a trend of 
increasing Be abundance with decreasing temperature (in Fig. 
[6]), as we would need to compensate for the weaker blend by 
increasing the Be abundance to fit the 3131 line. This increase in 
the Be abundance would also make the 3130 line much stronger 
than observed. Choosing the Fe blending line, we can get a good 
simultaneous fits of both the 3130 and 3131 lines. Therefore, we 
conclude here that the Fe line seems to be a better solution for 
the blend, although in the future it might still be proven that it is 
not the right one. 

Although the Be line is not the dominant line at low temper- 
atures, the blend should still be sensitive to the Be abundance, as 
long as the blending lines are well constrained. Although this 
is not perfectly true, from the exercise we conducted we be- 
lieve to have as good a control over the blends as is currently 
possible. Thus, using spectrum synthesis, Be abundances can be 
determined down to at least r e fj ~ 4700 K, which is the low- 
temperature end of our sample. This is illustrated in Fig. [5] where 
we show the spectra of our coolest sample star in comparison 
with synthetic spectra with different Be abundances. 

To finalize, two points must be stressed. First, for as long 
as the blends remain not well known, we can not claim to deter- 
mine absolute Be abundances at the low T e ff range with sufficient 
precision. The relative comparison between similar stars in the 
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Fig. 5. Fit to the region of the Be lines in the star IC 2391 VXR 
67 with r e jf = 4680 K. The solid line represents the observed 
spectrum, the dotted line a the best fit synthetic spectrum with 
A(Be) = 1.30, and the dashed lines represent spectra without Be 
and with + 0.25 dex from the best fit. 
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Table 4. Uncertainty of the Be abundance. 



□ IC 2391 
O IC 2602 



Fig. 6. Beryllium abundances as a function of the effective tem- 
perature of the stars in IC 2391 (open blue squares) and IC 2602 
(open red circles). The error bar in temperature is shown only on 
the upper left side of the figure. 



same T e g range, however, should be robust. Second, as the Fe 
line we adopt is the stronger among the possible choices, our Be 
abundances should be on the low-end. Therefore, our conclusion 
that there is no significant Be depletion in the stars analyzed here 
(see Section|4]below) is robust. 

3.4. Uncertainties in the Be abundances 

The Be abundances are affected by random uncertainties coming 
from the atmospheric parameters and by uncertainties in the de- 
termination of the pseudo-continuum during the spectrum syn- 
thesis. To estimate the effect of the atmospheric parameters, we 
can change each one by its own error, keeping the other ones 
with the original adopted values, and recalculate the abundances. 
Thus, we measure how the variation of one parameter affects the 
abundances. We assume the uncertainties introduced by each pa- 
rameter to be independent of the other ones. The uncertainty due 
to the continuum was determined by estimating the sensitivity of 
the derived Be abundance on the adopted continuum level. This 
uncertainty is mostly related to the S/N of the spectrum, and thus 
can be slightly different from star to star. Estimations of these ef- 
fects are listed in Table [4] The calculations were done for three 
stars, spanning the rage in temperatures of the sample (VXR67, 
R14, and SHJM2). In the plots that follow these uncertainties 
were considered as typical values also for the other sample stars 

with similar temperatures. 

As mentioned in Sec. |2l iRandich et al.l d2002l) has also an- 
alyzed one of sample stars, star SHJM2 of IC 2391. The at- 
mospheric parameters they used are the same adopted here. An 
abundance of A(Be) = 1.11 +0.13 was derived. Our own value, 
A(Be) = 1.20 + 0.13 agrees with theirs within lcr. 

3.5. Lithium abundances 

Lithium abundances for IC 2 391 and IC 2602 were fir st calcu- 
lated bv lStauffer etal] (Il989t) and IRandich et all (Il997l) . respec- 
tively. The samples of member st ars of these two clust ers were 
then reanalyzed and extended by IRandich et al.l d2001l) . The Li 
abundances listed in Tab. [3] were all taken from this last ref- 



Star 


°~ Tcff 


C loge 


°> 


O" [Fc/H] 




°~ total 


SHJM2 


±0.05 


± 0.10 


±0.00 


±0.03 


±0.05 


±0.13 


R14 


±0.05 


±0.10 


±0.05 


±0.00 


±0.10 


±0.16 


VXR67 


± 0.10 


±0.20 


±0.10 


±0.00 


±0.20 


±0.32 



Notes. The effect introduced by the uncertainties of each of the atmo- 
spheric parameters and by the uncertainty of the fitting itself related to 
the level of the continuum and the S/N. 
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Fig. 7. Ber yllium abundances as a function of the Li abundances 
derived by IRandich et al.l d2001h for the stars in IC 2391 (open 
squares) and IC 2602 (open circles). 



erence. Li abundances are particularly sensi tive to temperature 
and the T e s adopted in IRandich et al.l d2001l) are for a few stars 
slightly different from the ones adopted here. The differences 
amount to at most 70 K. According to the analysis presented in 
IRandich et al.ld2001l) . these differences could result in 0.10-0.15 
dex differences in the Li abundances. Differences that are within 
the uncertainties of the analysis itself (for our sample stars the 
uncertainties in A(Li) listed in IRandich et al.l d2001l) vary be- 
tween 0.18 and 0.30 dex). 

The Be abundances, on the other hand, are only weakly sen- 
sitive ^toJhe_T S £ ; r^causej of this, and because the Li abundances 
from IRandich et ail d2001l) are, anyway, consistent within them- 
selves, we decided to just adopt their values for A(Li) instead of 
recomputing the abundances. 



4. Discussion 

4.1. Mixing in the pre-main sequence 

Solar-type stars start their lifes as fully convective objects with 
large radii and low central temperature. As they contract, the 
stellar temperature rises and may eventually reach values high 
enough to burn the light elements. The rise in temperature can 
also lead to the development of a radiative core while the convec- 
tive layer retreats to the upper layers. The radiative core will ap- 
pear at later ages the smaller the mass and the higher the metal- 
licity of the star. Light element depletion will only be observed 
if the temperature in the convective zone reaches the values nec- 
essary for nuclear burning. If the temperature is high enough for 
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Fig. 8. Comparison with the Be abundances from selected results of the literature. In the left panel the abundances are shown as 
a function of T e s and in the right panel as a function of mass. As in the previous figures, our results for IC 2391, IC 2602, and 
the Hyades are shown as op en blue squares , open red circles, and green st a rred s ymbols, respectively. Abundances for additional 
Hyades stars collecte d from iRandich et aT] d2007h and iGarcia Lopez et al.l d 19951) are shown as full g reen triangles. Abundances 
for Hyades stars from lBoesgaard & King ( 20021) and from Pl eiades stars fromlBoesgaard et al.l d2003bl) are shown as "x" and "+" 
symbols, respectiv ely. Abundances for NGC 2516 stars from lRandich et aD d2007l) are shown as full black circles and abundances 
for field stars from lSantos et al.l d2004l) are shown as open triangles. The dashed line indicates the solar value. 



depletion only in the radiative core, no depletion is observed in 
the photosphere. 

Models and observations have shown that stars more mas- 
sive then ~ 1.20 M do not deplete Li during the PMS. The 
amount of depletion increases for less massive stars. At about 
0.6 M Q all Li will be burned before the radiative core is able 
to develop. As beryllium is less fragile than Li, it is expected 
to survive longer. Significant Be depletio n is expected by the 
models only for masses b elow ~ 0.6 M Q dBodenheimeri 1 1 966t 
IChabrier & Baraffall997l) . This expectation, however, has not 
been empirically tested before and this is our aim here. 

Before discussing the Be abundances, lets briefly recall the 
be havior of the L i abund ances in our two clusters. As discussed 
in IRandich et al.l d2001l) . stars more massive then ~ 1 .0 M© 
do not show signs of Li depletion, while cooler stars are de- 
pleted by an amount that increases with decreasing temperature. 
IRandich e t al. (2001) concluded that the late-G and early-K stars 
in IC 2602 present a scatter in Li abundances that is similar, but 
not as large, as the one seen in the Pleiades. IC 2391 did not 
show strong evidence of a scatter, but this could be caused by 
low number statistics. Some correlation with rotation is seen, 
but not as strong as in the Pleiades. 

The dependence of Li depletion on temperature was found to 
be similar in the two clusters down to ~ 5000 K. At lower tem- 
peratures, the stars in IC 2391 seem to be more depleted than the 
ones in IC 2602. Moreover, some of the stars in the two clusters 
are more depleted in Li than similar stars in the older Pleiades. 
These observations point that the pattern of depletion may be 
affected by the individual properties of each cluster. 

The observed correlation between Li and rotation in main- 
sequence stars in young clusters is a strong indication that the 
evolution of stellar angular momentum during the PMS affects 
the efficiency with which Li is burned. In addition to the increase 
in rotation due to the contraction of the star, the angular mo- 
mentum evolution is also affected by the interaction of magnetic 
fields with the disk and/or the stellar wind. 



Given the complexity of the processes involved, abundances 
of both Be and Li are useful to probe the range in temperatures 
achieved in the convective zone during the PMS. Thus, demon- 
strating whether the expectation from the models, i.e. that Be is 
not affected, is important. 

In Fig. [6] we show the Be abundances of the stars as function 
of the effective temperatures. The mean abundance of our stars 
is A(Be) = 1.23 + 0.07. No trend with temperature is detected. 
Taking into account that the spectra of many of our sample stars 
are affected by rotational broadening, the Be abundances we de- 
rive are expected to b e slightly overestimated (see discussion in 
ISmihanic et alj|20lol) . Therefore, given the uncertainties, there 
is no significant difference between the abundances measured in 
our stars and that observed in the Sun (where A(Be) = 1.10). 

Recently, Tak eda et al.l d201 lb proposed that most early-G 
stars suffer some degree of gradual Be depletion. This conclu- 
sion is based on a weak correlation found between A(Be) and 
v sin i in a large sample of solar analogues. However, the obser- 
vations of Be in old open clusters and the similarity of the Be 
abundances among cluster members of different ages and field 
stars seem to argue against that, suggesting instead that no Be 
depletion has taken place in stars between 6100 K < T e ff < 5600 
K. If that is indeed the case, then Fig.|6]strongly suggests that our 
cooler stars also did not suffer any Be depletion. Below 5000 K, 
we were only able to determine the Be abundance in one star, 
VXR67 (see Fig. [5]). Although its abundance has a larger error 
bar, it is fully consistent with the remaining objects. There is no 
indication of Be depletion. 

This conclusion is reinforced by Fig. [7] where we plot the 
Be abundances as a function of the Li abundances derived in 
IRandich etall d2001l) . In the case where Be is affected by PMS 
depletion, one would expect the abundances of Li and Be to cor- 
relate, i.e. that stronger Be depletion should be seen in stars with 
stronger Li depletion. As no such trend is seen in Fig.|7]there is 
again no indication of Be depletion. 
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We recall here that, as discussed in Sec. 13.31 the Fe line we 
adopt as a blend is the stronger among the possible choices. 
Therefore, any systematic effect introduced by this choice would 
result in a smaller Be abundance for cooler stars. Thus our con- 
clusion that there is no significant Be depletion in the stars ana- 
lyzed here is robust. 

4.2. Depletion on the main sequence 

In Fig. [8]we present a comparison of our results with Be abun- 
dances of other clusters and field stars. The references from 
where the abundances were compiled are given in the caption 
of the figure. All of them derived a solar abundance of beryllium 
similar to the one we obtain here, and thus are in a similar scale. 
Exception are the res ults of Boesgaard and collaborators who 
found A(Be) = 1.30 dBoesgaard et al.l2003bl) . If we were to ex- 
press the results using the Sun as reference, we would be in the 
same scale. The addition of their results in Fig.[8]do change our 
conclusion below, as their stars are located in the same tempera- 
ture range where all other results also do not find Be depletion. 

The masses of the stars shown in the left panel of Fig. [8] 
were interpolated from theoretical isochrones, for the age of the 
cluster and the given temperature of the star. The isochrones 
were calculated with the internet servefl dedicated to pre-main 
sequence and mai n-sequence model calculations presented by 
ISiess et al ] (120001) . From these isochrones, it seems that the 
masses of stars older than the ZAMS might be slightly overesti- 
mated (by -0.10 M ). 

The sequence from the clusters to the field in Fig. [8] can be 
seen as a sequence of age: from IC 2391 and IC 2602 (~ 50 Myr), 
to NGC 2516(~ 150 Myr) the Hyades (~ 600 Myr), to the field 
stars of ISantos etail d2004l) (older than ~ 1 Gyr). It is seen that 
for stars cooler than ~ 5600 K (mass lower than ~ 1.10 M Q ) 
the degree of Be depletion seems to increase with age. There is 
also a clear correlation between Be depletion and temperature 
(mass). The depletion is stronger the lower the mass. As there 
seems to be no Be depletion for stars at ages of 50 and 150 Myr, 
we conclude that the Be depletion in this temperature range is 
taking place d uring the main seq uenc e, in agreement with the 
conclusions of lSantos etafl d20041) and lRandich et al.1 d2007l) . 

4.3. Comparison with the models 

We now compare our abundances with expectations from the- 
oretical models. Isochrones for the e arly-main sequen ce were 
calculated using the internet server of ISiess et al.l d2000h for so- 
lar metallicity (Z = 0.02), with and without overshooting, and 
for 5 different ages (1, 10, 50, 125, and 600 Myr). In addition 
to the usual quantities (luminosity and temperature) the models 
offer additional information such as the surface abundance of the 
light elements (LiBeB). 

When looking at the plots comparing abundances and mod- 
els, it is important to remember that a star does not keep its 
temperature constant during the whole age interval for which 
we show models. In other words, a fixed temperature value in 
the following plots correspond to stars with different masses de- 
pending on the age. For example, at 1 Myr, a star with ~ 3 M Q 
has 5000 K, but at 50 Myr (onwards) a star with ~ 0.90 M Q has 
this temperature. 

In Fig. [9] we compare the Li abundances with the models 
without overshooting. Both IC 2391 and IC 2602 have ages of 
about 50 Myr while the Hyades has ~ 600 Myr. The plot shows 

7 http://www.astro.ulb.ac.be/~siess/server/iso.html 
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Fig. 9. C omparison of the Li abundances with the theoretical 
models of Siess et al. The curves show the Li abundance 

for stars with solar metallicity at the main sequence as a func- 
tion of effective temperatures. Five different ages are shown: 1 
and 10 Myr as black solid curves, 50 Myr as a red dotted curve, 
125 Myr as a green short-dashed curve, and 600 Myr as a blue 
long-dashed curve. In addition to the stars of IC 2391, blue open 
squares, and IC 2602, red open circles (IC 2602), we show two 
stars from the Hyades as green starred symbols. 
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Fig. 10. Comparison of the Be abundances with the theoretical 
models of Siess et al. 1 J2000I) . The curves show the Be abundance 
for stars with solar metallicity at the main sequence as a function 
of effective temperatures. Four different ages are shown: 10 Myr 
as a black solid curve, 50 Myr as a red dotted curve, 125 Myr as 
a green short-dashed curve, and 600 Myr as a blue long-dashed 
curve. Symbols for the stars are as in Fig. [9] 



what was discussed before, i. e., the models predict too much 
Li depletion for cool stars and no depletion during the main 
sequence. The difference between the young clusters and the 
Hyades, however, suggest that some depletion took place dur- 
ing the main sequence in this older cluster. 

The same comparison is done with the Be abundances in Fig. 
[TOl The models show slightly higher Be abundance because of 
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a difference in the choice of reference solar abundance. Apart 
from that, the observations agree very well with the models. In 
the temperature range investigated, which for an age of 50 Myr 
correspond to a mass range between 0.80 < M/M G < 1 .20, no Be 
depletion is expected and none is seen. 

Formally, the inclusion of the Hyades stars in these plots and 
their comparison with the Z = 0.02 models could be questioned 
as they are more metal rich ([Fe/H] = +0.13). From the mod- 
els point of view, metallicity is ex pected to influence the de- 
gree of light eleme nts depletion (e.g. lPiau & Turck-Chiezd2002t 
ISestito et al.ll2.QQ6h . An increased opacity would affect the con- 
ditions at the bottom of the convective zone during PMS. In 
iRandich etail d2007l) . however, models of PMS Be depletion 
were calculated for both the solar and the Hyades composition 
for a star of 0.9 M B . It was shown that the effect caused by 
metallicity in this case is negligible. Thus, the observed differ- 
ence between the Be ab undances can not b e explained by the 
difference in metallicity. ISestito et al.l (120061) has indeed inves- 
tigated in detail the effect of heavy-element opacity on PMS Li 
depletion. They concluded that changes in the abundances are 
not enough to result in agreement between models and Li ob- 
servations in low-mass stars of young clusters. Observationally, 
no significant difference between Li depletion patterns has been 
found among clusters with similar ages and different metallic - 
ities, argu ing against a strong metallicity depe ndence of PMS 
depletion (I Jeffries & Jameslll999tlJeffriesll2006l) . 

Similar comparisons are shown in Figs.[TTIand[T2lbut against 
the models including overshooting. The discrepancy with the Li 
data is shown to increase when overshooting is taken into ac- 
count. For Be however, the inclusion of overshooting is still not 
enough to cause Be depletion in the range of parameters of our 
sample. From these plots we see that we were still not able to 
reach down to effective temperatures where Be depletion during 
the pre-main sequence phases could be detected. If one can push 
the determination of Be abundances to even lower temperatures, 
it would be possible to look for the effect of overshooting on Be 
depletion. Such measurement a re, however, be very c hallenging. 

As discussed previously in IRandich et al.l (120071) . the small 
difference seen between the Be abundances in the Hyades and in 
IC 2391/2602, seems to point to yet another shortcoming of the 
models. The mixing mechanism depleting Li during the main 
sequence seems also to be able to deplete Be, although not as 
dramatically. 

4.4. Evolution of Be in the disk 

The average Be abundance in the two clusters is A(Be) = 1 .24 
± 0.09 and 1.22 + 0.06, for IC 2391 and IC2602 respectively. 
These values can be compared to the average of other young 
clusters to search for variations of the Be abundance in the disk 
and to investigate signs of Be evolution with time in the Galaxy. 
For this comparison we selected from the literature only stars 
with 5400 K < r eff < 5800 K, i.e. stars that show no signs of Be 
depletion. For the clusters analyzed by Boesgaard and collabora- 
tors the averages are: A(Be) = 1.33 ± 0.06 for 10 stars of the 600 
Myr old Hyades; A(Be) = 1 .26 + 0.07 for 3 stars of the 120 Myr 
old Pleiades; A(Be) = 1 .28 for 2 stars of the 500 Myr old Coma; 
A(Be) = 1. 25 for 1 stars of the 300 Myr old Ursa Major mo v- 
ing g roup dBoesgaard & Kingl 120021: iBoesgaard et alJ l2Q03allbL 
2004). All these average values seem to be in agreement with 
each other, given the uncertainties. Although one should keep in 
mind that there is a small spread in metallicity between the clus- 
ters (from [Fe/H] = -0.09 for Coma and [Fe/H] = +0.13 for the 
Hyades) and that the results are not exactly on the same scale as 
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Fig. 11. C omparison of the Li abundances with the theoretical 
models of Siess et al. including overshooting. Lines and 

symbols are as in Fig. [9] 



ours, when seen together they indicate a lack of strong variations 
in the initial Be abundance among the clusters. 

By comparing the Be abundance in these young clusters with 
the solar abundance, one can search for evidence of Be enrich- 
ment in the disk. This task if made difficult by a possible dis- 
agreement between the solar photospheric and meteoritic abun- 
dances, A(Be) = 1.10 and A( Be) = 1 .42 respective l y. Thi s mete- 
oritic value recommended by lAnders & Grevessel (Il989l) seems 
to indicate a photospheric depletion by a factor of two. 

iBalachandran & Belli dl998l) and lBell et all (1200 ll) . however, 
argued that by accounting for a possible missing near-UV con- 
tinuum opacity the photospheric Be abundance increases to 
A(Be) = 1.41, and thus is brought into agreement with the me- 
teoritic one. This conclusion is corroborated by a re-analysis of 
the Be photospheric so lar abundance using 3D hydrodynamical 
models (lAsplundl2004 . wher e a recommende d va lue of A(Be) = 
1.38 is found. As discussed in lAsplundl (120051) and lAsplund et al.l 
d2009h . Be abundances in the Sun seem to be largely insensitive 
to NLTE and 3D effects and this increase is only due to the miss- 
ing opacity. It is interesting however to notice that two recent re- 
determinations of the Be abund ance in a CI chondrite meteorite 
yields a value of A(Be ) = 1.32 dMakishima & Nakamurall2006l 
see also lLoddersll20 1 Ol) . which is smaller than the currently rec- 
ommended photospheric value. The implications of this smaller 
value have not been discussed in the literature yet, but this lower 
value might suggest that the 'missing opacity' has been overesti- 
mated in the above analyses. In their an alysis of Be abundances 
in solar analogues, iTakeda et al.l d201 ll) reached similar conclu- 
sion. 

As discussed before, the Be abundances we derived in the 
young clusters show no evidence of PMS depletion and are in 
agreement with the one we derived for the Sun. Thus we can 
draw the two following conclusions: (1) there is no evidence 
for significant Be depletion in the Sun, in agreement with the 
conclusions of IRandich et al.1 (f2002l l2007t) : (2) the similarity of 
the abundances also indicate that there has been no significant 
Galactic enrichment of Be in the last ~ 4.5 Gyr since the Sun 
was formed. 
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Fig. 12. C omparison of the Be abundances with the theoretical 
models of Siess et al. including overshooting. Lines and 

symbols are as in Fig.fTOl 

5. Summary 

Beryllium abundances were calculated for ten confirmed mem- 
bers of the young pre-main sequence clusters IC 2391 and IC 
2602. We also recalculated Be abundances for two main se- 
quence stars of the Hyades. With a careful abundance analy- 
sis, we were able to provide actual detections for fast-rotating 
stars and for stars with effective temperatures down to ~ 4700 
K. These young stars (~ 50 Myr) have just arrived in the main 
sequence. Therefore, any alteration in the surface abundance of 
the light elements Li and Be should have taken place during the 
PMS phase. 

We show that all the stars analyzed have, within the uncer- 
tainties, the same Be abundance even though they have different 
degrees of Li depletion. This result confirms empirically for the 
first time the theoretical prediction that Be abundances are not 
affected by mixing during the PMS for stars in the range between 
0.80 < M/Mq < 1.20, as expected by the models. In addition, as 
discussed in the literature, for stars less massive than the Sun the 
mixing mechanism depleting Li during the main sequence seems 
also to be able to deplete Be. 

At even lower temperatures (~ 4000 K), models with and 
without overshooting differ in the predicted amount of Be de- 
pletion. It is not clear, however, if the Be abundances can be 
properly determined in such cool stars. The attempt to push the 
determination of Be abundances to even lower temperatures, al- 
though very challenging, could provide important information 
to differentiate between these models, testing the importance of 
overshooting during the PMS. 
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